Eigenvalue correlations and the distribution of ground state angular momenta for 

random many-body quantum systems 



J. Barea 

Center for Theoretical Physics, Sloane Physics Laboratory, Yale University, 
P.O. Box 208210, New Haven, Connecticut 06520-8120, U.S.A. * 



o 
o 

(N 

< 



o 

=5 



(N 
> 
O 



oo 
o 



R. Bijker and A. Frank 

Instituto de Ciencias Nucleares, Universidad Nacional Autonoma de Mexico, A. P. 70-543, 04510 Mexico, D.F., Mexico 

(Dated: April 16, 2009) 

The observed preponderance of ground states with angular momentum L = in many-body 
quantum systems with random two-body interactions is analyzed in terms of correlation coefficients 
(covariances) among different eigenstates. It is shown that the geometric analysis of Chau et al. 
can be interpreted in terms of correlations (covariances) between energy eigenvalues thus providing 
an entirely statistical explanation of the distribution of ground state angular momenta of randomly 
interacting quantum systems which, in principle, is valid for both fermionic and bosonic systems. 
The method is illustrated for the interacting boson model. 
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I. INTRODUCTION 

Low-lying spectra of many-body quantum systems of- 
ten display a high degree of order and regularity. In the 
case of atomic nuclei, despite their complexity and the 
large number of degrees of freedom involved, they of- 
ten exhibit simple features, such as pairing properties, 
surface vibrations and rotational motion in even-even 
nuclei. Conventional wisdom is that regularities arise 
from symmetries of the Hamiltonian, which lead to in- 
variances that severely constrain the many-body motion. 
While some of these symmetries are exact (e.g. rota- 
tional and reflection invariance), others are approximate 
(e.g. isospin). These global symmetries, however, do not 
explain by themselves the regular patterns observed. Fur- 
ther assumptions about the nature of the nucleon-nucleon 
interaction are required. Thus, a strongly attractive pair- 
ing force between like nucleons has been shown to be re- 
sponsible for the remarkable constancy of the excitation 
energy of the first excited L = 2 states in the Sn isotopes, 
while deformation and rotational behavior is known to 
arise from an attractive quadrupole-quadrupole interac- 
tion between neutrons and protons These striking 
patterns as well as many other correlations have been 
shown to be robust features of low-energy nuclear behav- 
ior, which signal the emergence of order and collectivity. 
In every case the patterns arise as a consequence of par- 
ticular forms of the nucleon-nucleon interaction. Most 
features of low-lying nuclear spectra have thus been ex- 
plained in terms of a short-range pairing interaction and 
a long-range quadrupole force. 

Hence, it came as a big surprise when Johnson, Bertsch 
and Dean [2] found that for even-even nuclei shell-model 
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Hamiltonians with random two-body interactions (the so- 
called Two-Body Random Ensemble or TBRE) are very 
likely to yield a ground state with angular momentum 
L = 0. In fact, the probability for a L = state to 
become the ground state turns out to be much larger 
than that expected on the basis of the fraction of L = 
states in the model space. Similar regularities were later 
found to exist as well in bosonic [3| and electronic Q 
many-body quantum systems with random interactions, 
so the occurrence of these regular phenomenona in spite 
of the random nature of the two-body interactions seems 
to be a rather generic feature. The unexpected results of 
Johnson et al. are reminiscent of other statistical results 
associated with correlated sets of numbers, which are per- 
ceived as contrary to expectations. This is the case, for 
example, for Benford's Law, which deals with the coun- 
terintuitive relative frequency distributions of digits in a 
given data set and which is related to scale invariance [H]. 

The observation of a statistical preference of L = 
ground states for ensembles of random two-body interac- 
tions has sparked a large number of investigations to fur- 
ther explore the properties of these random systems and 
to understand the mechanism for the emergence of reg- 
ular ordered spectral features from random interactions 
0, 0, H, [§] . The appearance of ordered spectra in systems 
with chaotic dynamics is a robust property, that does not 
depend on the specific choice of the (two-body) ensem- 
ble of random interactions 0, [l(J EH time-reversal 
symmetry (lot ], the restriction of the Hamiltonian to one- 
and two-body interactions [T3|, nor is it limited to yrast 
states with small angular momentum L = 0,2,4 [14| as 
used in the original studies 0, 0]. Despite the progress 
and new insights in understanding the appearance of or- 
dered spectra from random interactions, the proposed ex- 
planations offer partial solutions of the problem, without 
being able to account for all observed regular phenomena 
in systems of randomly interacting fermions and bosons. 
Among others, we mention induced pairing [l5j], geomet- 
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ric chaoticity for randomly interacting fermions 0, EB], 
mean-field analysis for the interacting boson model and 
the vibron model [T3|, spectral widths [13, [HI and an 
empirical method based on the eigenvalues of each inde- 
pendent two-body matrix element fl2| . 

The robustness of the numerical results for both 
systems of randomly interacting fermions and bosons 
strongly suggests that an explanation of the origin of the 
observed regular features has to be sought in the many- 
body dynamics of the model space and/or the general 
statistical properties of random interactions, a conclu- 
sion which is suggested by many theoretical studies. In 
this respect, the approach of Chau et al. for diagonal 
Hamiltonians stands out, since it makes use of the geom- 
etry of the model space, is valid for both fermions and 
bosons, and allows to calculate the ground state proba- 
bilities exactly. 

In this article, we show that the geometric analysis 
of [19] can be interpretated in terms of correlations (co- 
variances) between energy eigenvalues to provide an en- 
tirely statistical explanation of the distribution of ground 
state angular momenta of randomly interacting many- 
body quantum systems. In addition, we show that the 
method can be extended to non-diagonal Hamiltonians 
by using perturbation theory. 



II. COVARIANCES 

In [iji], spectroscopic properties of quantum systems 
with random interactions were given a geometric inter- 
pretation. In particular, it was shown that diagonal 
Hamiltonians, i.e. whose energy eigenvalues depend lin- 
early on the two-body matrix elements, can be associated 
with a geometric shape (convex polyhedron) defined in 
terms of coefficients of fractional parentage and/or gen- 
eralized coupling coefficients. The probability for a cer- 
tain state to become the ground state is then related to 
the angles at the vertices. In this approach, geometry 
arises as a consequence of strong correlations implicit in 
many-body quantum systems. Random tests can be un- 
derstood in this context as sampling experiments on this 
geometry. 

Let us consider many-body quantum systems for which 
the energy eigenvalues e, can be expressed as linear com- 
binations of two-body matrix elements r m (as is the case 
for any interaction between identical fermions with j < | 
or identical bosons with I < 3) 

d 

m— 1 

The coefficients c l m contain the information on the many- 
body quantum system via angular momentum coupling 
coefficients, coefficients of fractional parentage, etc. The 
statistical properties of such a system can be studied by 
taking an ensemble of random two-body matrix elements 



r m in which the random variables r m are chosen indepen- 
dently on a Gaussian distribution with zero mean and 
unit width. The covariance coefficients between two en- 
ergies are given by 

i N 

M - j^YlW -WW- <&)) 

d 

— ^ , c m c m = c ' c3 ' J (2) 
m— 1 

where N is the size of the ensemble and p indicates the 
p-th realization of the ensemble. The covariances depend 
on the relative angle between the vectors c 1 and c- 7 . 

It was shown by Chau et al., that all energies are con- 
fined to a convex polytope (i.e. a convex polyhedron in 
d dimensions), and that only the states that are located 
at the vertices of this polyhedron can become the ground 
state [iji]. The probability for a state at vertex j to be- 
come the ground state depends on the angle @jfi 
where the sum is over all faces that contain the vertex j 
and 9jf is the angle subtended at vertex j in the face /. 
For d = 2 and d = 3 dimensions the explicit forms are 

j 2 2tt ' 

From these considerations it follows directly that the 
probability that a certain state can become the ground 
state is related to the covariances of the energies. Sup- 
pose that the angle 9jf at vertex j in the face / is de- 
fined by the vertices i, j and k of the convex polygon. A 
straightforward application of the cosine rule shows that 

(ej) + (e l e k ) - (e l e j ) - (eje k ) 
COS Vjf = — = — — 
^(e?) + <e?) - 2(e i e j ) y /(el) + (e|) - 2(e 3 e k ) 

(4) 

We thus find that the probability for a state which is 
located at a vertex j of the convex polytope to become 
the ground state is determined by the covariance coeffi- 
cients of the energy eigenvalues. It is important to note 
that this relation is exact and may be used to explain 
the app roximate results based on the spectral widths 
[ToL 1 1 Sl ] - Eq. (j4]) is valid for any many-body quantum 
system whose energies are linear functions of the random 
variables and holds for both bosons and fermions. 



III. RESULTS 

The relation between the distribution of ground state 
angular momenta and the covariances between energy 
eigenvalues of Eq. Q provides a completely statistical 
interpretation of the distribution of ground state angu- 
lar momenta of randomly interacting many-body quan- 
tum systems. To the best of our knowledge, this is the 
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first time such a connection has been derived in explicit 
form. Since the geometric method of Chau et al. is valid 
for diagonal interactions for which the energy eigenval- 
ues depend linearly on the random interactions, Eq. Q 
provides a statistical interpretation of the exact results 

of [H. 

For systems in which one has both diagonal and off- 
diagonal matrix elements the method of covariance co- 
efficients cannot be applied directly, since the energy 
eigenvalues show a nonlinear dependence on the random 
variables. In this section, we show that also for non- 
diagonal systems the distribution of ground state angu- 
lar momenta can be obtained by treating the off-diagonal 
interactions in perturbation theory. In this way, one ob- 
tains an approximate expression of the energy eigenval- 
ues which is linear in the random variables and hence is 
amenable to a statistical analysis based on the covari- 
ances. Obviously, this is only valid as far as the distri- 
bution of ground state angular momenta is concerned. 
Realistic ground state wave functions have a far more 
complicated structure that goes way beyond the use per- 
turbation theory. 

As an illustration, we discuss two schematic Hamilto- 
nians of the Interacting Boson Model (IBM), a model for 
collective excitations in medium and heavy mass nuclei 

As a first example, we consider the IBM Hamiltonian 



cos 



which describes a transition between vibrational (sin x — 
0) and 7-unstable nuclei (cosx = 0). In the general case, 
the Hamiltonian has to diagonalized numerically to ob- 
tain its eigenvalues. The eigenstates can be labeled by 
the total number of bosons N, the boson seniority v and 
the angular momentum I. The ground state has either 
v = 0, v = 1 or v = N. In the present study the angle x 
is taken on the interval — ir < x < 7T, so that all possible 
combinations of attractive and repulsive interactions are 
covered. For this reason, the Hamiltonian effectively de- 
pends on two random coefficients cosx an d sin% which 
can be either attractive or repulsive. 

The distribution of ground state angular momenta was 
obtained exactly in a mean-field study fT3] - For even 
values of the number of bosons N, the ground state has 
v = in 75 % of the cases and v — N in the remaining 
25 %. For odd values of N, the ground state is either 
v = 0, v = 1 or v = N for 50 %, 25 % and 25 % of the 
cases, respectively. 

The Hamiltonian of Eq. ([5]) has a spherical ground 
state for sin y < | cosxl and a deformed one for sinx > 
I cos xl [13, EH- In the first case, the appropriate basis 
is that of the (7(5) limit of the IBM, \N,n,v). The first 
term in the Hamiltonian of Eq. J5]) is diagonal in this ba- 
sis, whereas the second term contains both a diagonal and 
a nondiagonal contribution. The latter interaction can be 
treated in first order perturbation theory so that the ap- 
proximate energy spectrum has a linear dependence on 



Table I: Distribution of ground state probabilities P v (in %) 
for N — 15 and N = 16 bosons, obtained exactly (Full) and 
in the diagonal approximation (Diag) using Eq. for the 
spherical solution and ([7]) for the deformed case. 







Spherical 


Deformed 


Total 


N 
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Full Diag 


Full Diag 


Full Diag 


15 





50.0 48.6 


50.0 50.0 


50.0 49.0 




1 


33.3 34.7 


0.0 0.0 


25.0 26.0 




N 


16.7 16.7 


50.0 50.0 


25.0 25.0 


16 





83.3 83.3 


50.0 50.0 


75.0 75.0 




1 


0.0 0.0 


0.0 0.0 


0.0 0.0 




N 


16.7 16.7 


50.0 50.0 


25.0 25.0 



the parameters 



E nv = cos x rid + 



■[(N-n){N-n- 1) 



smx 
A(N- 1)' 

+(n - v)(n + v + 3)] . (6) 



For the deformed solution, the appropriate basis is that 
of the 50(6) limit of the IBM, \N,a,v). In this case, the 
second term of Eq. <[5j) is diagonal, whereas the first term 
contains diagonal and off-diagonal contributions. In first 
order perturbation theory, the energy spectrum is given 

by 



E av = sinx 



(N - a)(N + a + 4) 
4(AT- 1) 



where Q a = (N,n,v\N, a, v) 0. 

In Table HI we show a comparison of the exact results 
and the ones obtained in first order perturbation the- 
ory, using the energy eigenvalues in Eqs. © and (0, 
and the correspondence between the covariances and dis- 
tribution of ground state probabilities of Eq. |(4]). The 
results obtained in perturbation theory are remarkably 
close to the exact ones. There are some slight deviations 
for the spherical solutions, whereas in the deformed re- 
gion there is an exact correspondence for the distribution 
of ground state angular momenta between the diagonal 
approximation and the full calculation. 

The Hamiltonian of Eq. ((5j) is a special choice, since in 
addition to the number of bosons and the angular mo- 
mentum, also the boson seniority v is a conserved quan- 
tum number. Let us now consider a Hamiltonian in which 
this is not the case. As a second example, we take a 
schematic IBM Hamiltonian that describes the transi- 
tion between vibrational and rotational nuclei. In the 
notation of [23(], this Hamiltonian is given by 

H 2 = cosxd t -d+^Y[(2s t s t -d t -d t )(h.c.) 

+(2sW±T/7dW)W ■ (h.c.)j . (8) 

The Hamiltonian of Eq. ((Sj) has a spherical minimum 
for — 7T < x < arctan i and a deformed minimum for 
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Figure 1: (Color online) Percentage of ground states with 
angular momentum L = (circles), and L = 2 (triangles) 
for 10 < N < 16 bosons interacting via the Hamiltonian of 
Eq. with — ty < \ < n calculated exactly (solid lines) and 
in first order perturbation theory (dotted lines). 



arctan | < x < 7r. For x = arctan |, the system exhibits 
a first-order phase transition between spherical and de- 
formed nuclei [2ll. |23|. Just as in the previous example, 
the distribution of ground state angular momenta can 
be obtained exactly in a mean-field study. For N = 3k 
(k = 1,2,.. .) the ground state has L = in 75 % of the 
cases and L = 2N (the maximum value of the angular 
momentum) in the remaining 25 %. For all other values 
of the total number of bosons N, the ground state has 
either L = 0,L = 2orL = 2N for 50 %, 25 % and 25 % 
of the cases, respectively. 

In order to test the correspondence between the co- 
variances and the distribution of ground state angular 
momenta, we treat the Hamiltonian of Eq. JH]) in first or- 
der perturbation theory, in which the first term is treated 
as a perturbation to the second in the deformed region 
and vice versa in the spherical region. In Fig.HJ we show 
a comparison of the percentages of ground states with 
L = and L = 2 obtained exactly and in first order per- 
turbation theory. Also in this case, the agreement is very 
good, although not at the level of precision as shown in 
TableUfor the Hamiltonian of Eq. © . The deviations ob- 
served for the probabilities for ground states with L = 
and L = 2 may be reduced in second order perturbation 
theory. The percentage of ground states with the maxi- 
mum value of the angular momentum L = 2N is 25 % in 
all cases, both exactly and in perturbation theory. 

In both examples, the distribution of the ground state 
angular momenta is reproduced to a great level of preci- 
sion with the present procedure, which is a combination 
of perturbation theory to get energy eigenvalues which 
are linear in the random two-body matrix elements, and 
the relation that we derived between the covariances and 
the probabilities of the ground state angular momenta. 



The energies of many-body quantum systems are 
strongly correlated for the simple reason that they are 
all eigenvalues of the one and the same Hamiltonian. 
These strong correlations may modify the distribution 
of ground states from naive expectations, as was the case 
for the preponderance of L = ground states observed 
for calculations in the nuclear shell model and the IBM 
with random interactions 0. Due to the strong corre- 
lations, it should not be surprising that the probability 
distribution of the ground state angular momenta for an 
ensemble of random two-body interactions is related to 
the covariances between energy eigenvalues. 



IV. SUMMARY AND CONCLUSIONS 

In this manuscript, we addressed the problem of the 
emergence of regular features in many-body quantum 
systems with random interactions. As an example, we 
focussed on the probability distribution of ground state 
angular momenta in nuclear models, especially the un- 
expected observed dominance of L = ground states 
@, S|- Even though many authors have tackled the sub- 
ject, no general explanation for the preponderance of 
L = ground states has been found. A general feature of 
these systems is that the energy eigenvalues are strongly 
correlated, since all many-body matrix elements are ex- 
pressed in terms of a relatively small number of random 
two-body matrix elements. 

In this article, we proposed to explore the correlations 
between energy eigenvalues (or covariances) in more de- 
tail. We established, to the best of our knowledge for 
the first time, an explicit relation between the proba- 
bility distribution of ground state angular momenta and 
the covariances between energy eigenvalues. For diago- 
nal (but in no way trivial) Hamiltonians, our formulas 
in terms of covariance coefficients are exact. For non- 
diagonal Hamiltonians it is not possible to give closed 
expressions but, by means of perturbation theory, nearly 
exact results can be found for the ground-state percent- 
ages. 

This new relation provides a purely statistical inter- 
pretation of the abundances of ground states with L = 
observed in numerical studies of many-body systems with 
random interactions which is, in contrast to all other ex- 
planations, valid for both fermions and bosons (e.g. the 
nuclear shell model and the IBM in nuclear physics). As 
an example of this procedure we studied two schematic 
nondiagonal Hamiltonians of the IBM and found an ex- 
cellent agreement between the approximate results ob- 
tained in perturbation theory and the exact ones. Al- 
though these results were illustrated in the context of a 
specific model, we believe them to be sufficiently gen- 
eral to propose an entirely statistical explanation which 
is valid in general for any many-body quantum system. 
In future work, we will apply the present method to prob- 
lems in the nuclear shell model. 

As a final comment, an explanation for the emergence 
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